Long ranged interactions in carbon atomic chains 
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Based on first-principles calculations we revealed fundamental properties of infinite and finite size 
monatomic chains of carbon atoms in equilibrium and under an applied strain. Distributions of bond 
lengths and magnetic moments at atomic sites exhibit interesting even-odd disparity depending on 
the number of carbon atoms in the chain and on the type of saturation of carbon atoms at both 
ends. It was shown that, the 7r-bands of carbon atomic chains behave as a one dimensional free 
electron system. A local perturbation created by a small displacement of the single carbon atom at 
the center of a long chain induces oscillations of atomic forces and charge density, which are carried 
to long distances over the chain. These long ranged oscillations are identified as Friedel oscillations 
showing 1/r decay rate in one dimensional systems. 

PACS numbers: 63.22.-m, 73.21.Hb, 73.20.Mf 



I. INTRODUCTION 

Carbon atomic chains are one-dimensional (ID) al- 
lotropic form of carbon atom, which has also allotropic 
forms in different dimensionalities, such as 3D diamond 
and graphite, 2D graphene, quasi ID nanotube and 
quasi 0D fullerens. sp^-hybrid orbitals are indigenous to 
the dimensionality (D=l,2,3) of these allotropic forms. 
The tetrahedrally coordinated sp 3 -bonding stabilizes the 
open diamond structure. The sp 2 -bonding together with 
7r-bonding maintain the planar stability of honeycomb 
structure of graphene and attributes several exceptional 
physical properties.^ Covalent bonding of sp D=1 hybrid 
orbitals along the chain axis together with 7r-bonding 
of perpendicular p x and p y orbitals are responsible for 
the linear stability of the chain. Earlier carbon atomic 
chains and their functionalized forms (CACs) have been 
investigated intensively, despite the lack of consensus on 
whether they can really be synthesized PHU These stud- 
ies have predicted a wide range of interesting properties, 
which can make CACs a p oten tial material for future 
nanotechnology applications!^ 

Linear CAC structures have either identical double 
bonds, called cumulene or alternating short "single" and 
strong "triple" bonds, called polyyne as described in 
Fig. [IJa). While cumulene is metallic with a quan- 
tum ballistic conductance of Ae 2 /h due to two degen- 
erate, half-filled p x and p y bands crossing the Fermi 
level, it is vulnerable to Peierls instability.^ Hence, 
through the displacement of alternating carbon atoms by 
5 = 0.018 A the unitcell is doubled and a band gap of 0.32 
eV at the edge of Brillouin zone (BZ) is opened to lower 
the total energy per atom by 2 meV. While segments 
of polyynes terminated with H atoms at both ends (H- 
C rl -H), have been produced^ up to considerable lengths 
(n=20), cumulene production is relatively difficult due 
to their frailty. Small cumulene chains ter minate d by 
H 2 , (H 2 -C n -H 2 ) groups have been synthesized! 10 1 11 1 Free- 
standing CACs from graphene flakes were produced^ 
by using energetic electron irradiation inside a transmis- 



sion electron microscope (TEM). Concomitantly, it was 
demonstrated that not only CACs, but also SiACs and 
BNACs can be derived from their corresponding honey- 
comb structure under uniaxial tensile stress in the plastic 
deformation ranged Much recently, polyyne structure 
consisting of 44 carbon atoms have been produced.^ 

In this work, we predict that the even-odd dispari- 
ties are attained in the distributions of bond lengths and 
atomic magnetic moments of finite size CACs depend- 
ing on the type of saturation of carbon atoms at both 
ends. Even more remarkable is that a local perturbation 
through the displacement of a single chain atom creates 
atomic force and charge density oscillations, which prop- 
agate to long distances in the chain. We show that, these 
long ranged couplings in CACs can be explained in terms 
of Friedel theoryPSl 



II. METHODS 

Our predictions are obtained from the state-of-the- 
art first-principles plane wave calculations carried out 
within the density functional theory (DFT) using 
PAW potentials.^ The exchange correlation potential 
is approximated by generalized gradient approximation 
(GGA) using PBE functional. A plane-wave basis set 
with kinetic energy cutoff of 450 eV is used. All CACs 
are treated by supercell geometry and atomic positions 
and lattice constants are optimized by using the conju- 
gate gradient method, where the total energy and atomic 
forces are minimized.^ The vacuum separation between 
the CACs in the adjacent unit cells is taken to be at least 
10 A. The convergence for energy is chosen as 10~ 5 eV be- 
tween two steps, and the maximum Hcllmann-Feynman 
forces acting on each atom is less than 0.01 eV/A upon 
the ionic relaxation. Phonon dispersions were obtainned 
using the force constant method with forces calculated in 
a (40 x 1 x 1) supercellPl 
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FIG. 1. (Color Online) (a) Schematic representations of or- 
bitals and structure of cumulene (left) and polyyene (right), 
(b) The energy per unitcell of two carbon atoms versus the 
dimerization A = S/c c , that is defined as the ratio of displace- 
ment of one of the carbon atoms from the cumulene positions, 
S — di — d,2, to the lattice constant of cumulene, c c . In each 
curve the minimum energy is set to zero, (c) Variations of 
band gap, E g , and tension force, Ft, with applied strain, e. 
(d) Calculated phonon dispersions of polyyne structure. 



III. INFINITE CARBON ATOMIC CHAINS 

The variation of energy as a function of dimerization A 
in Fig. [ljb) , shows that cumulene structure corresponds 
to a metastable state which transforms to polyyne struc- 
ture without facing any energy barrier. This instabil- 
ity is also reflected in the phonon modes of cumulene 
structure. The longitudinal acoustic modes of cumulene 
attain imaginary frequencies near the zone boundary. 
The analogy of this transition can be found in graphene, 
where 1.5 fold bonds of the equilateral hexagons break 
into single and double bonds. It was shown that, if the 
non-local part of the exchange-correlation is taken into 
account through hybrid functionals, the high symmerry 
phase of graphene becomes considerably less stableP2"23 
Thus, we have performed a structural optimization of cu- 
mulene and polyyne structures also using HSE, B3LYP 
and PBE0 hybrid functionals, as they are implemented 
in VASPp2H231 Indeed, the energy difference between the 
cumulene and polyyne structures, which is 2 meV ac- 
cording to PBE functional, is increased to 47 meV, 80 
meV and 104 meV when HSE, B3LYP and PBE0 hybrid 
functionals are used, respectively. Interestingly, while the 
ground state of infinite CAC is polyyne and hence cumu- 
lene transforms to polyyne, cumulene structure, by itself, 
can be stabilized by charging. We deduced interesting 
effects of charging on infinite and finite CACs, which is 
beyond the scope of the present study. 

In Fig. [ljb) one can see that applying stress to the 
CACs enhances dimerization both energetically and spa- 



tially. As the lattice constant is increased by 8% the en- 
ergy difference between polyyne and cumulene structures 
is increased to 40 meV as compared to 2 meV difference 
in the absence of stress. Also the spatial dimerization 
relative to the lattice constant is 3 times higher in 8% 
strained structure than that in the unstrained one. The 
increase of dimerization with strain results in the increase 
of the band gap, as shown in Fig. [l|c). The band gap 
reaches its maximum value of 2.87 eV (which is an order 
of magnitude higher than that in the unstrained case) at 
a critical strain e c =0.18. At this critical strain, the band 
gap switches from direct to indirect, since the minimum 
of the second conduction band starts to dip below the 
first conduction band. This rapid change of band gap 
with strain is interesting and can make CACs a potential 
candidate for strain gauge type nanodevice applications. 

Fig. [ljc) also presents the variation in the tension 
force Ft with strain. The tension force is defined as 
Ft = —8Et /dep, where Ft is the total energy per unit 
cell. Ft curve starts with a linear region at small strain 
values, which corresponds to the elastic regime. Follow- 
ing this region, the slope of this curve, which is propor- 
tional to the sound velocity, starts to decrease.^ This 
lasts until the sound velocity drops to zero where the 
tension curve reaches its maximum. This is the critical 
point occurring at e c , beyond which CAC cannot sustain 
the long wavelength perturbations. 

The phonon dispersions of unstrained polyyne is pre- 
sented in Fig. [ijd). Degenerate transversal phonon 
branch of polyyne structure mimics the folded version of 
that of cumulene, except a small gap which separates this 
branch into acoustical and optical modes. High frequency 
(wave number) longitudinal branch of polyyne also re- 
mind the folded version of that of cumulene but now 
there is a dramatic difference in a sense that imaginary 
frequencies at long wavelength disappear. The transver- 
sal acoustic mode loses its quadratic behavior near the 
r point under strain. Longitudinal acoustic and optical 
modes, however, change dramatically with strain. The 
dispersion of the optical branch decreases with increasing 
strain. In the limit of very high strains where individual 
carbon dimers do not interact, this mode is expected to 
converge to a flat line with frequency corresponding to 
the vibrational motion of an isolated carbon dimcr. The 
sound velocity related with the slope of the longitudi- 
nal acoustic mode decreases with increasing strain. The 
imaginary frequencies appear near the T-point beyond 
e e , which indicates the onset of instability under the long 
wavelength perturbations. 

The phonon modes are derived by using the direct 
method where the force constants are calculated in a 
(40 x 1 x 1) unitcell comprising 80 atomsP^Such a long 
unitcell comprising 80 atoms was necessary because force 
constants in longitudinal directions decay rather slowly. 
On the other hand, the transversal force constants decay 
rapidly, which leads to appearance of quadratic terms 
in the phonon dispersions of transversal modes near the 
center of BZ. It is interesting to compare the decay of 
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FIG. 2. (Color Online) The even-odd disparity: Distribution of C-C bond lengths and magnetic moments, fii, at chain atoms 
of a finite CAC including n C atoms, which is saturated by m hydrogen at both ends, namely (H m -C„-H m ). m=0 corresponds 
to the bare CAC with free ends. H-saturated (m=l,2,3) CACs with even n and bare CACs with odd n have nonmagnetic 
ground state. All other CACs have total magnetic moment of 2/j.b distributed on the atomic cites as shown in the right panel. 
The dotted lines are a guide to the equilibrium bond distances of cumulene and polyyne structures. 



force constants of CACs with those of graphene. In the 
honeycomb structure the magnitude of the out of plane 
(transversal) and in-plane force constants of 5 th neigh- 
bor, which is 4.3 A apart, is about 1/100 and 1/ 50 of 
the nearest neighbor force constants, respectively! 26 * 27 ! 
Similar to graphene, in CACs, the transversal force con- 
stant is reduced to ~ 1/100 at the 4 th neighbor which is 
5.1 A apart. However, the longitudinal force constant of 
CACs is reduced to about 1/50 only at the 22 nd neighbor, 
which is 28.2 A apart. This slow decay in the magnitudes 
of longitudinal force constants of CACs is accompanied 
with sign oscillations starting from 3 rd neighbor. The 
slow decay of longitudinal force constants in CACs indi- 
cate long ranged nature of specific interactions, which is 
the principal subject of our study and will be treated in 
the forthcoming parts. 



IV. LONG RANGED INTERACTIONS IN 
FINITE CARBON ATOMIC CHAINS 

Infinite CACs cannot exist in reality, but the proper- 
ties of long finite CACs are expected to converge to those 
of infinite ones. There are two new parameters which can 
affect the properties of finite CACs dramatically. These 
are the number of carbon atoms n forming the chain and 



the end effects. To simulate the end effects we either 
leave both ends of CACs bare or passivate each of them 
by one, two or three hydrogen atoms. A similar analy- 
sis, which confirms our results, was performed earlier for 
only bare CACs and those passivated with one hydrogen 
from both edgesP^ We also expect that the CAC passi- 
vated with two (three) hydrogen atoms from both ends 
to behave similarly to CACs attached between graphene 
sheets (diamond blocks) Figure [2] presents the distri- 
bution of bond lengths and atomic magnetic moments 
of finite CACs having n=60 and 61 carbon atoms satu- 
rated by m hydrogen atoms from both ends. It turns out 
that for a given m the structural and magnetic proper- 
ties can be classified in two classes depending on whether 
n is even or odd. This means that the distribution of 
bond length and magnetic moments of (H2-C60-H2) and 
(H-C61-H) structures, for example, is similar to that of 
(H2-C100-H2) and (H-C101-H), respectively. These trends 
leading to even-odd disparity are confirmed for chains 
with n = 20, 21, 60, 61, 100, 101 and for m = 0, 1, 2, 3. 

When edges are passivated by only one H atom, the 
C atoms at both edges make single bond with these H 
atoms. This forces the type of C-C bond to the adja- 
cent chain atom to be triple. Then the next C-C bond is 
forced to be single and so on. This alternation of triple 
and single bonding is also reflected to the short and long 
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FIG. 3. (Color Online) Forces and difference charge density perturbations formed in finite Cioo, B50N50 and AI100 atomic chains 
due to a small longitudinal displacement, 6z, of the 50 th atom near the center of the chain. The upper panel is a schematic 
representation of the geometry of the system. The middle panel presents force distribution for three nearest neighbors (left) and 
extensions up to 35 th neighbor. In the similar manner the bottom panel presents the difference charge density perturbations. 
Atomic magnetic moments (which are not shown in the figure) also behave similarly, and have long ranged oscillations in CAC. 



bond length alternations. A situation contrary to this oc- 
curs, when the edges are passivated with three H atoms. 
This time H atoms are arranged tetrahedrally and the 
outermost C-C bonds at both edges is forced to be single 
so that the bond length alternation starts with a longer 
bond. These structures are so similar that the C-C bonds 
in the structure passivated with one H atom remain un- 
changed if these H atoms passivating both ends are re- 
placed by CH3 groups. The bond length alternations in 
these structures are enhanced at the edges. For even 
values of n, the middle parts acquire a bond length al- 
ternation between two values corresponding to the bond 
lengths of polyyne. For odd values of n, however, the 
bond length alternation originated from both edges com- 
pete at the middle parts and the pattern presented in 
Fig. [2] is formed. 

When the ends are passivated by two H atoms, the type 
of C-C bonds at the ends becomes double. This time all 
other C-C bonds acquire the same type. As a result, for 
sufficiently long CACs, the bond length alternation in the 
middle parts become negligible while the bond lengths 
are very close to that found in cumulene. Carbon atomic 
chains making two bonds with cone-terminated carbon 
nanotubes from both edges were shown to have cumulene 
type bonding, which corroborates our findings .^Follow- 
ing the similarity observed between the CACs passivated 



by single and triple H atoms, the bare CACs and CACs 
passivated by two H atoms have similar atomic structure. 
That is, bare CACs also have double bonds between C 
atoms and have negligible bond length alternations in the 
middle parts. It is remarkable that an atom at the center 
of a chain as long as 100 C atoms, is affected by the type 
of passivation at its edges. 

Similar trends are obtained in the distribution of 
atomic magnetic moments, where their values oscillate 
until long distances. Since the effect of a specific type 
of saturation occurring at both ends of a long chain can 
be carried over to distant neighbors, resulting even-odd 
disparities imply a long ranged nature of couplings in 
CACs. 

We now consider the aforementioned atomic force 
and difference charge density oscillations of bare Cioo, 
B50N50, B 10 o and Al 10 o chains generated by a slight dis- 
placement of the 50th atom. Despite its zigzag metal- 
lic ground statep^ here we consider linear Al-chain for 
the sake of comparison. After the structural relaxation, 
all atomic positions are kept fixed except the 50 th atom, 
which was displaced in longitudinal direction by Sz taken 
to be 0.02 A for Cioo, B50N50, Bi 00 and 0.04 A for Alioo 
chain. Difference charge density is obtained by subtract- 
ing the charge density of the perturbed structure from 
that of the unperturbed one, which is averaged in planes 
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perpendicular to the chain axis. In the left panel of Fig.[3j 
the force and linear charge density perturbations of only 
three nearest neighbors of the 50 th atom is shown. The 
right panel of the Fig. [3] presents the extensions of these 
perturbations up to the 35 th neighbor of the 50 th atom. 
One can see that these extensions are negligible and show 
no obvious pattern in Al chains. Similar pattern is also 
seen in Boron atomic chains. In CACs, however, the per- 
turbation extensions are long ranged and exhibit a decay- 
ing oscillatory behavior. The envelope of these oscillatory 
extensions fits to a 1/r decay rate (r being the distance 
to the 50th atom) for both force and linear charge den- 
sity perturbations. Similar pattern is also seen in B50N50 
chain, but both force and charge density oscillations are 
weaker compared to that of Cioo- 



V. DISCUSSIONS AND CONCLUSIONS 

The above oscillatory decay in the linear charge den- 
sity perturbation provokes an association with so called 
Friedel oscillations.^ Inserting an impurity charge to 
an electron gas results in the accumulation of electronic 
charge, which screens the Coulomb potential introduced 
by the impurity. Since the wave vector of the electron 
gas is limited by its density, the accumulated charge ac- 
quires decaying Friedel oscillations. For sufficiently large 
separations, r, from the impurity, the Friedel oscillations 
are proportional to the function sin(2kpr) /r D , where kp 
and D stand for the Fermi wave vector and dimension- 
ality of the electron gas system, respectively.^] Recently 
it was shown that, upon inclusion of an impurity in a 
two-dimensional graphene structure, where D = 2, the 
charge density oscillations proportional to sin(2k F r)/r 2 
are obtainedP^In case of CAC, one can think of the dis- 
placement of the 50 th atom near the center of the chain 
as the insertion of two impurity charges having opposite 
signsPS Since the carbon chain system is one dimensional, 
the Friedel theory can explain the 1/r decay rateP3 In 
cumulene structure the double degenerate 7r-bands are 
half filled, so the Fermi wave vector is kp = it/2c c . In 
polyyne structure bands are folded and the Brillouin zone 
is halved and the Fermi wave vector of filled 7r-bands is 
kp = ir/cp. Thus the periodicity of observed charge den- 
sity oscillations fits the Friedel theory for both polyyne 
and cumulene structures (and also for BN chains). 

One should note that, Friedel theory holds for free elec- 
tron like systems. In our case, the 7r-bands of atomic 
chains are expected to behave as free electron systems, 
since their charge density have nodes on the axis where 
the ions lie. Figure [4] presents the band structure of A1A1, 
BB, BN and polyyne chain structures. All structures are 
considered in a unitcell having two atoms for the sake of 
comparison. The left panel presents the whole band pro- 
file, while the right panel zooms to the 7r-bands of chain 
structures and compare them with free electron disper- 
sion. As seen in the right panel of Fig. |4j the 7r-states 
of all chain structures considered here have a dispersion 
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FIG. 4. (Color online) Energy band structures of A1A1, BB, 
BN and polyyne chain structures. A closer view of 7r-bands 
together with the dispersion profile of one-dimensional free 
electron gas system is presented in the right panel. Zero of 
energy is set to the minima of 7r-bands for each structure. For 
the sake of clarity, the bands of the chain structures beyond 
the first Brillouin zone in the extended scheme are not shown. 



profile which is very similar to that of a free electron. In 
fact the effective masses of 7r-bands of BB, CC, A1A1 and 
BN chain structures exceed the mass of free electron by 
only 2%, 7%, 8% and 17%, respectively. Here the effec- 
tive masses are calculated from the curvature of the E(k) 
profile around the T point. 

BB and A1A1 chains both have six valence electrons 
in a unitcell consisting of two atoms. The folded dou- 
ble degenerate 7r-bands of these structures are half filled, 
which result in kp = ir/Ac. Considering that, the 7r-band 
dispersion of BB chain is closer to free electron disper- 
sion compared to that of polyyne, one expects the BB 
and A1A1 chains to have Friedel oscillations with a peri- 
odicity of 4c. However, there is no such profile present 
in force or charge density oscillations of A1A1 chains, as 
seen 111 Fig. [3] We think that, in A1A1 and BB chains 
the Friedel oscillations, which is a feature of a free elec- 
tron systems, is suppressed by ionic potential. In CC 
and BN chains the oscillations have a periodicity of 2c, 
which resultes in very small charge density perturbations 
near the atomic sites. In A1A1 and BB chains, however, 
the charge density oscillations with a period of 4c are ex- 
pected to have maxima on atomic sites. This enhances 
the contribution of ionic potential, so that, 7r-bands de- 
viate from free electron like behavior. As a result, the 
charge density perturbations, being a feature of free elec- 
tron like system, are suppressed and the system finds a 
self consistent energy minimum with no Friedel oscilla- 
tions. Note that, force and charge density oscillations are 
weaker in BN chains compared to that of CACs. This is 
because, the 7r-bands of BN chains considerably deviate 
from free electron behavior, as seen in Fig. [4j 

The force distribution observed in Fig.|3]can be related 
to the charge density oscillations; as mentioned before, in 
the absence of any perturbation, the minimum of energy 
is met by symmetric charge densities around the atomic 
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sites. As the charges start to move in one direction, the 
atomic sites start to feel a force in the same direction. To 
the first order, these two events are directly proportional 
and this explains the oscillating 1 jr decay of the atomic 
forces. Note that, the 50 th atom and its neighbors feel 
the elastic forces generated by displacement. Thus, they 
are not in the same direction as the charge accumulation. 

In conclusion, a local displacement of an atom in car- 
bon atomic chains creates long ranged oscillations of 



atomic forces and charge density. These remarkable and 
fundamental features of carbon chains are explained in 
terms of Friedel Theory developed for one-dimensional 
free electron system. 
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